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The CAST 'GPS 4000 simulator can create scenarios using external trajectories. The 
following information must be provided: time from a given epoch (a constant sampling rate must 
be used); position, velocity and acceleration in Earth Centered Earth Fixed (ECEF) coordinates; 
elements of the matrix that transforms from the body-axis coordinate system to the ECEF frame; 
and the angular velocity of the body-axis system relative to ECEF coordinates. The initial 
latitude, longitude, altitude, and UTC time must also be provided during the scenario setup. 

The simulator recomputes the positions and velocities using the given accelerations and a 
constant jerk model. The results are the internal “truth file”. The basic algorithm provided by 
CAST for each ECEF component is 
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The jerk over each interval is computed by simple differencing. 
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Substituting this into the basic algorithm gives 
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The external trajectory files require the linear accelerations, but this data is not available 
from radar or GPS systems and must be derived from the position or velocity data. The 
technique used at Kennedy Space Center to obtain the accelerations is very simple and is based 
on the philosophy that the measured data should not be changed unless absolutely necessary, for 



example, if it is very noisy or has gaps. Noisy data is smoothed with a simple moving average 
filter whose width is chosen based on the rapidity of the dynamic changes and the data rate, but 
usually about ±10 data points is sufficient. Any gaps are filled in with quadratic interpolation. 
If only position data is available, the velocity is obtained from the position by midpoint 
differencing, 
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The acceleration is obtained from the velocity in a similar manner, 
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Although this technique has worked for us, it must be used carefully because the 
derivatives can be very noisy unless the original data is smooth and well-behaved. There are 
many other ways of smoothing the data and obtaining the derivatives that may be more desirable 
or appropriate for a well-informed user. 

Of course, one wants to know how well the simulator reproduces the position and 
velocity of the original input file. To answer this question, point-by-point comparisons were 
made between the positions and velocities from the internal truth and external trajectory files for 
three operations: the Athena Kodiak Star launch from Kodiak, Alaska; a fictitious three-stage 
modified Minuteman missile launch from Wallops Island, Virginia called the Wallops Express; 
and an actual weather balloon flight from Wallops Island, Virginia. The results are shown on the 
following pages. The velocity and acceleration profiles are also given. All times are seconds 
from launch. 

The two launch trajectories contain less than 10 minutes of data, and the balloon flight 
contains about 1 hour and 45 minutes of data. The velocities in the external trajectory files were 
converted to an East North Up system to agree with the output format of the internal truth file. 

No attitude data was available and only the center-of-mass coordinates were used. 

The Kodiak and Wallops Express external trajectory files were based on theoretical 
nominals with clean, smooth position and velocity data that was quadratically interpolated from 
10 Hz to 50 Hz before numerically differentiating the velocities as described above. Because the 
data was very smooth, this higher rate resulted in smaller jerks. The Wallops balloon external 
trajectory file was based on position-only 1 Hz data collected from an on-board GPS receiver 
during the actual flight. This data contained noise and random movements of the balloon and 
was quadratically interpolated to 1 0 Hz before numerically differentiating the position and 
velocities: interpolating to a higher data rate resulted in unacceptably large jerks. No smoothing 
was done on any of the data. 

The rocket scenarios required no ramp-up of the velocities, but the balloon scenario 
required a two-second ramp of the initial velocities because there were large jerks when the 
balloon was released. 

The position differences between the input and internal truth data are less than a few 
meters and the velocity differences are typically much less than a few centimeters per second. 

The largest differences occur during times of large or rapidly changing accelerations (i.e., jerks). 
Because numerical integration is used, the mean errors gradually increase. For all but the most 
demanding applications, it is reasonable to conclude that that these differences are negligible. 
Although is not known why the balloon data is noisy at about T+4000 s, atmospheric turbulence 
is a likely cause. 
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